Pull counting statistics of Andreev scattering in an asymmetric chaotic cavity 
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We study the charge transport statistics in coherent two-terminal double junctions within the 
framework of the circuit theory of mesoscopic transport. We obtain the general solution of the 
circuit-theory matrix equations for the Green's function of a chaotic cavity between arbitrary con- 
tacts. As an example we discuss the full counting statistics and the first three cumulants for an 
open asymmetric cavity between a superconductor and a normal-metal lead at temperatures and 
voltages below the superconducting gap. The third cumulant shows a characteristic sign change as 
a function of the asymmetry of the two quantum point contacts, which is related to the properties 
of the Andreev reflection eigenvalue distribution. 

PACS numbers: 74.50.-|-r,72.70.+m,73.23.-b,05.40.-a 
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I. INTRODUCTION 

Mesoscopic heterojunctions with superconductor (S) 
and normal-metal or semiconductor (N) leads exhibit a 
rich variety of phenomena that have been studied through 
conductance and noise for some time already^^^'^ ; see also 
Refs. 14-18. In the absence of disorder charge transport 
through an interface between S and N can be described 
by the Blonder-Tinkham-Klapwijk (BTK) model/ while 
the effects of scattering by impurities and interface rough- 
ness can be taken into account using the scattering the- 
ory of mesoscopic transport. Recent experimental 
advancements"'^^"^^ raised interest in the statistical prop- 
erties of the charge transfer in mesoscopic systems. Cur- 
rent noise contains information on the temporal corre- 
lation of quasiparticles relevant for transport and orig- 
inates from the fluctuations of quasiparticle occupation 
numbers or random scattering at barriers or impurities. 
Further, it can be used to probe internal energy scales 
of the system or the effective charge of the elementary 
transport mechanisms, "'^^"^^ as well as to detect the cor- 
relations intrinsic to the many-body state of entangled 
systems. -'^^ 

If the charge transfer events are uncorrelated, the zero- 
frequency noise power at low temperatures (shot noise) 
has the Poisson value Pi = 2el, where e is electron 
charge and / is the average current through the sam- 
ple. This is the case for the charge transport in a low- 
transparency normal-state tunnel junction at zero tem- 
perature and low bias. In less-opaque normal-state junc- 
tions shot noise is usually suppressed below the Pois- 
son value, with the deviations characterized by the Fano 
factor F = P//2e/ = E„?^«(l " ?^«)]/E„r„, where 
Tn are transmission eigenvalues of the junction. For a 
metallic diffusive wire^^^^" and for an open symmetric 
chaotic cavity^^ Fano factors have universal values^^'^^"'^" 
F = 1/3 and F = 1/4, respectively, which do not depend 
on microscopic properties like the impurity concentration 
or geometry of the sample. 

The effect of superconductivity leads to a doubled shot 
noise in subgap transport through low-transparency S/N 
tunnel junctions^** and in diffusive normal wires in con- 



tact with a superconductor. ^^'^■^ This doubling can be 
understood as a consequence of the effective charge dou- 
bling in the Andreev process. '^^ In general, however, the 
noise is also affected by the change of transmission prop- 
erties of the structure due to the superconductor prox- 
imity effect and the doubling is not generic. For exam- 
ple, in an open symmetric cavity with a superconducting 
and normal-metal lead the Fano factor has the valuc"'^^ 
F « 0.60 which is more than two times larger than in the 
corresponding normal-state junction (F = 1/4), in agree- 
ment with recent experimental results. At bias voltages 
on the order of the superconducting gap both normal and 
Andreev scattering processes contribute to transport'^^''^^ 
and the picture of the effective charge carriers fails. '^^ 

The statistical theory of mesoscopic transport, "^^ full 
counting statistics, provides the complete characteriza- 
tion of the charge transfer and has led to a new and fun- 
damental understanding of quantum transport phenom- 
ena in nanoscale conductors. •^^■'^^ Higher-order moments 
of the charge transport statistics may provide additional 
information to the current noise measurements. ^^'^^''^^ 
The third-order correlations of voltage fluctuations across 
the nonsuperconducting tunnel junctions have been mea- 
sured by Reulet et al?^ and recently by Bomze et al.,'^'^ 
the latter confirming the Poisson statistics of the charge 
transfer at negligible coupling of the system to environ- 
ment. The semiclassical theory of higher-order cumu- 
lants based on the Boltzmann-Langevin equations has 
been developed recently by Nagaev et al.^^ The quantum- 
mechanical approach to full counting statistics based 
on the extended Keldysh-Green's function technique'*°'''^ 
in the discretized form of the circuit theory"*^'^^ was 
put forward for multiterniinal circuits by Nazarov and 
Bagrets.^"* This approach can describe junctions with dif- 
ferent connectors and leads as well as multiterminal cir- 
cuits in unified and very general way. Within the circuit 
theory, the doubling of shot noise (i.e., F — 2/3) was 
found in diffusive S/N junctions in both the fully coher- 
ent eV ^ Eth and the completely incoherent eV ^ Exh 
regimes,''^''*® where Etii is the inverse diffusion time. At 
intermediate bias voltages eV Exh , the transport is af- 
fected by electron-hole coherence leading to an enhance- 
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ment of the differential Fano factor. ^^■''^ The influence of 
coherence effects on the full counting statistics and noise 
in other single and double S/N junctions was studied in 
Refs. 48-50. The current correlations in a three-terminal 
chaotic cavity operated as superconducting beam splitter 
were studied by Samuelsson and Biittiker,^^ yielding un- 
usual positive cross correlations. Prom a study of the full 
counting statistics these positive correlations can be at- 
tributed to the uncorrelated injection of Cooper pairs. 

In this article we study the full counting statistics of 
coherent charge transport in a chaotic cavity using the 
circuit theory of mesoscopic transport. We show that the 
system of matrix equations for the Green's function of 
the cavity can be solved, effectively, as a system of scalar 
equations independently of the type of leads and without 
resorting to the matrix components or parametrizations. 
As an application we find the Green's function for an 
open asymmetric cavity between arbitrary leads. For the 
special case of a cavity between the superconductor and 
normal metal, we find the cumulant generating function 
and the first three cumulants and discuss the interplay 
between superconducting proximity effect and scattering 
properties of the junction. The results are compared with 
those for a normal-state junction^^ and for different cou- 
plings of a cavity to the leads. Current correlations in 
a structure with high-quality contacts between a cavity 
and superconductor have been studied experimentally by 
Choi et al."^^ recently. 



II. MODEL 

The system we consider is a chaotic cavity coupled to 
two leads by mesoscopic junctions characterized by the 
transmission eigenvalues {T^^} and {T,!^''}, respectively. 
Charging of the cavity is negligible if the cavity is large 
enough and the conductances of the junctions §1^2 are 
much larger than the conductance quantum qq = jh. 
We assume an isotropic quasiparticle distribution func- 
tion inside the cavity due to chaotic scattering. The deco- 
herence effects as well as the energy dependence of trans- 
mission eigenvalues can be neglected if the total dwell 
time in the cavity is small with respect to time scales set 
by the inverse temperature and bias voltage. We apply 
the circuit theory of mesoscopic transport with the spe- 
cific parts of the system represented by the corresponding 
discrete circuit elements, as shown in Fig. 1. The leads 
are characterized by known quasiclassical matrix Green's 
functions Gi,2 which depend on the quasiparticle energy, 
lead temperature, chemical potential, and counting field 
and satisfy the normalization condition G\ = G\ = \. 
The lead Green's functions are not necessarily the ones 
of bulk electrodes - for example, they can be nodes that 
are part of a larger circuit. In the following we will re- 
fer to these nodes as leads. The formulation below is 
independent of the concrete matrix structure, provided 
the 'lead' Green's functions obey the normalization con- 
ditions. The chaotic cavity is represented as an internal 



node associated with an unknown Green's function Gc, 
which will be obtained from the matrix current conser- 
vation and the normalization condition G\ = \. Left and 
right junctions (i = 1,2), depicted as connectors, carry 
matrix currents^^ 



' h ^ 



'^Tn^ \Gi, Gc 



4 + Ti^^({G„GJ-2)' 



(1) 



which flow from the cavity into the leads. The current 
conservation Ii + I2 = for the Green's function Gc of 
the cavity reduces to 



[piGi+p2G2,Gc] =0, 



with 



rp{i) 

n 



^ 4 + Ti'\{Gi,Ge}-2) 



(2) 



(3) 



Here we have used the commutation of Pi{2) with Gi(2) 
and Gc, which is a consequence of the normalization 
G? = Gc = 1, and the matrix property 



A' 



1 



[A,{AB}]=o. 



(4) 



We can solve Eq. (2) while assuming that pi and p2 
depend only on the anticommutator of the Green's func- 
tions of the leads, pi = pi{{Gi, G2}) (see the Appendix), 
and commute with Gi, G2, and Gc in accordance with 
Eq. (4). As a result, the Green's function of the cav- 
ity can be expressed in terms of Green's functions of the 
leads in the form 



Gc = ^Gi + ^G2, 
c c 



(5) 



where the matrix c = c({Gi, G2}) accounts for the nor- 
malization of Gc- From Eq. (5) and by using the nor- 
malization conditions Gp = Gi = G| = 1 we obtain the 
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FIG. 1: A chaotic cavity coupled to the leads by two junctions 
with transmission eigenvalues {T,^^} and {t/i^'} (top). The 
discrete circuit-theory representation of tfic system is shown 
in the lower part. The leads and cavity are characterized by 
the corresponding matrix Green's functions. The junctions, 
depicted as connectors, carry conserved matrix currents. 
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system of equations 



C = Pi + P2 + PiP^y , 

002 = 2P2 + PlG , 



(6) 

(7) 
(8) 



where Gi = {Gi,Gc} and Q = {Gi,G2}- This system 
can be treated, effectively, as a system of scalar equations 
because all matrices that appear in Eqs. (6)-(8) depend 
only on Q and commute with each other. 

The cumulant generating function S{x) of charge 
transfer can be obtained as a sum of the actions of the 
connected pairs of nodes*^: 

^W=-| EE/ Trln(l + ^(ft-2)), (9) 

i=l,2 n '' 

where the total measurement time to is much larger than 
the characteristic time scale on which the current fluc- 
tuations are correlated. The x-independent term in the 
cumulant generating function [given by the normaliza- 
tion requirement S{x = 0) = 0] is omitted for brevity 
throughout this article. Also, we recall that because of 
current conservation, it is sufiicient to introduce a count- 
ing field X at one lead only, while the full counting field 
dependence can be obtained by setting X = Xi ~ X2- 

The cumulant generating function S depends only on 
the anticommutator Q (x) and is invariant to the exchange 
of the leads Gi ^ G2 or, equivalcntly, to exchange of the 
junctions {T^^^} ^ {T^^'} [see Eqs. (2) and (3)]. There- 
fore, the same invariance persists"'^'"'*'^^ in all coherent 
(or low bias) transport properties of two-terminal double 
junctions - such as current (conductance), noise (Fano 
factor), and higher cumulants - independently of the type 
of leads or specific properties of the junctions. This in- 
variance does not hold in the presence of dephasing,^*'^^ 
which can be modeled by an additional lead that carries 
the coherence leakage current. ^^"^'^'^^ 

In the following we consider an analytically tractable 
case of a chaotic cavity coupled to the leads by two quan- 
tum point contacts with A''i and N2 open channels, re- 
spectively. The transmission eigenvalues of the contacts 



are T^^ = 1 for n = 1, . . . , iVj and T, 
Prom Eq. (3) and Eqs. (6)-(8) we find 



otherwise. 



c = 
and 

Pl,2 
C 




4:Ni_N2 



(iVi + iV2)2 g + 2 



(10) 



1 / iVi-iVz 

2 I N1+N2 



4:NiN2 g -2 



(iVi + N2)^ g + 2 



The Green's function of the cavity and cumulant gener- 
ating function are given by Eqs. (5) and (9), respectively. 



A formally similar result has been obtained recently by 
Bulashenko''* using 2x2 Green's functions in Keldysh 
space which can be expanded over the Pauli matrices. 
Physically, this implies that the whole circuit is in the 
normal state, although it is permitted that the leads are 
nodes of a larger (normal-state) mesoscopic network. In 
our approach we do not rely on this expansion and make 
no assumptions on the particular matrix structure, ex- 
cept for the usual normalization condition. Therefore, 
Eqs. (10) and (11) are valid for any type of leads. For ex- 
ample, one lead can be superconducting, with the Green's 
function having Keldysh-Nambu matrix structure, or the 
chaotic cavity can be a part of the larger multiterminal 
network which consists of different heterojunctions. Ad- 
ditional degrees of freedom - for instance, spin - can be 
included as well. We emphasize, again, that our solution 
only resorts to the normalization condition of the leads. 
In the case in which the cavity is part of a larger network. 
Green's functions Gi_2 have to be determined by circuit 
rules. The result for Gc is valid in this case as well, which 
can simplify the numerical solution of larger circuits. 

It is interesting to check that an alternative approach 
can give the same result. Coherent connectors in the 
circuit theory are described by a cumulant generating 
function of the form 



— JdsJ^' dTp{T) Trln (l + ^[^(x) 



2] 



(12) 



where p{T) is the distribution of transmission eigenval- 
ues {r„} for the composite junction. Using the distri- 
bution of transmission eigenvalues for an open chaotic 
cavity,i6.42 



Pc(T) 




T-To 



(l-T)(l-ro; 



(13) 



for To < T < 1 and Pc(T') = otherwise, with Tq = 
(A^i - N2f/{Ni + N2f, we obtain S{x) as given by 
Eqs. (9) (11). This demonstrates the consistency of the 
circuit-theory approach with the random matrix theory 
of scattering matrices. 



III. SUPERCONDUCTOR - CHAOTIC-CAVITY 
- NORMAL-METAL JUNCTION 

In the following we calculate the statistics of charge 
transport through a chaotic cavity sandwiched between 
a superconductor and a normal-metal lead. We present 
a detailed analysis of the first three cumulants - cur- 
rent, current noise power, and the third cumulant of 
the current - for an open chaotic cavity at temperatures 
and bias voltages well below the superconducting gap A, 
when Andreev scattering is the dominant process of the 
charge transfer. At low energies and temperatures the 
4x4 matrix Green's function of the superconductor is 
Gs = Gi = 1 ai in the Keldysh(")(g)Nambu(~) space. The 
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Green's function = Gn{£,x) = G2 of the normal- 
metal lead incorporates the counting field according to 

Gjv = e-*4^^ G^e'^^^, (14) 

where fx = fi&s and fj and ai are Pauli matrices. The 
bare Green's function of the normal-metal lead is given 
by 



-CT3 



K 



1-2/ 







1-2/ 



(15) 



where /^ = {c:x.p[{±£+eV)/kBT] + l}^^ accounts for the 
voltage bias of the normal-metal lead, with the energy £ 
measured in respect to the chemical potential of the su- 
perconductor. From Eqs. (9)-(ll), we find the following 
expression for the cumulant generating function: 

Six) = -'-^ f d£ J2 
■' 3=1,2 

xln[r, + -^r2-64iV/(l + a)], (16) 



where 



1(2) 



a(7Vi-7V2)' + (3A^i(2)+A^2(i)) 



+ -^(1 + a)[a(7Vi - N2Y + (iV2 + iV| + 6iViiV2)2] . 

(17) 

Here, a = (e^'x _ 1)/+/- + (g-^'x _ 1)/+/- and /± = 
1 — /^, with a being related to the double-degenerate 
eigenvalues oi Q = {Gs,Gjv} given by \i.2 = ±.2is/a. 
From Eq. (16) we obtain the average current, the cur- 
rent noise power, and the third cumulant according to 
/ = i{e/ta)d^S\^=o, Pi - (2eVio)9^^lx=o, and Ci = 
~i{e^ /t^i)^^S\^=o^ respectively. They 



are 



/at fNfN)' 



(18) 



Pi =2Gs J d£ [if+fj, + /+/^) - 71 (/^/^ - Z^/^)'], 

(19) 

and 



Ci = 2eGs 



J d£ (/+ 



/at fNfN) 



X [1 - (7;^/,^ + f^ffN) + '^l2CfVfN - fi^fN?' 



with the total conductance of S-cavity-N junction, 

= + iV.)(l-^^ 

h \ q 

and the low-temperature Fano factor^^ 
l&NlN^iNi + N2) 



Fs 



(20) 



(21) 



(22) 



Here q ^ V(N\TN2fTWJ^, 71 = 1 - Fs/2, and 
72 = l-Fs[l~2NiN2{2Ni+N2)iNi+2N2)/q'^]. After 
the energy integration in Eq. (18), the usual relation / = 
GsV is obtained, while the integration in Eqs. (19) and 
(20) yields 



Pr 



= l + (Fs/2) [vcot\i{v)-l] 

^{l + {Fs/&)v\ \v\^l, 
\l-Fs/2 + {Fs/2)\v\, |t;|>l. 



(23) 



(24) 



and 



Cr 



eGsksT 



=12(71 - 72) 



+ v 



coth(ti) 
1 - 72 



coth^(w) 



(25) 



.3(71 - 72) 

[±12(71 - 72) 4(1- 371 272) w, ±w>l, 

(26) 

respectively, with v = eV/kBT. 

The cumulant generating function. Eq. (16), takes into 
account the superconductor proximity effect in the qua- 
siclassical approximation as well as interchannel mixing 
inside the cavity in the presence of quantum point con- 
tacts. In comparison with the normal-state junction (see 
Ref. 38), Eqs. (18) (20) contain products of electron and 
hole distribution functions due to the Andreev process, 
which is the mechanism of the charge transport. For ex- 
ample, the term f^fj^ = f^i^ ^ ./w) be interpreted 
as the probability for an electron emerging from the lead 
N to be reflected back as a hole, where and are 
the electron- and hole- state occupation numbers. The 
energy-independent prefactors Gg, 71, and 72 are also 
modified by the electron-hole correlations introduced by 
the superconductor. This change of transport properties 
due to superconductor proximity effect can be revealed 
by considering a general S/N junction with transmission 
distribution p(T), at low temperatures and bias voltages 
(ksT, \eV\ « A). In this case Eq. (12) reduces to 



Six) 



dRA PAiRA)Hl + RAaix)]. 



(27) 

where Ai_2 = ±2i^/a are double-degenerate eigenvalues 
oi Q. In the fully coherent regime, which we consider 
here, the distribution of Andreev reflection eigenvalues 
PAiRA) is simply related to the distribution of trans- 
missions p(T) of the corresponding normal-state junc- 
tion: PAiRA) = 2piT)dT /dRA. The probability of the 
Andreev reflection is given by Ra = T^/i2 — T)^ and 
appears in Eq. (27) due to electron-hole symmetry at 
energies well below the superconducting gap and the in- 
verse dwell time in the cavity. Normal scattering pro- 
cesses arc suppressed, and Six) depends on the count- 
ing field through e^"^^^, which accounts for the effective 
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charge of 2e that is transferred across the structure in 
each elementary event of Andrcev scattering. In the 
case of strong electron-liole dephasing within the struc- 
ture, electrons and holes decouple and the system can 
be mapped onto an effective normal-state junction, 
for which the cumulant generating function is given by 
Eq. (12) with the corresponding modification of trans- 
mission distribution p{T) and boundary conditions Q. 
In the crossover regime, transport through the structure 
is not simply related to the normal-state transmission 
properties and can be described by an effective energy- 
dependent distribution*^ pA which takes into account the 
effects of dephasing (coherence leakage currents) at char- 
acteristic energies on the order of Thoulcss energy. 

Expanding Eq. (27) in the x field and taking the 
derivatives, we obtain the current, current noise power, 
and the third cumulant in the coherent regime as given 
by Eqs. (18)-(20), with conductance 

Gs = 2 [ dTp{T) Ra, (28) 
Jo 

71 = 2Gs' dTp{T)Rl, and 72 = 2Gs' dTp{T)R% 
where Gs = Gs/(2e^ /h). In particular, the Fano factor 
and the slope of the third cumulant at high bias are given 
by 

4 /-i 

Fs = ^ dTp{T) Ra{1 - Ra) (29) 
Gs Jo 

and 

= ^ t dTp{T) Ra{1 - Ra){1 - 2Ra), (30) 
o{e 1) Us Jo 

respectively. These expressions are similar to the normal- 
state ones except for the effective charge doubling and the 
presence of the Andreev reflection probability Ra instead 
of normal transmission T, in agreement with previous 
results obtained within the scattering approach.^ Using 
the transmission distribution pdT) given by Eq. (13), we 
recover the results for an open asymmetric cavity, which 
were obtained from the circuit theory without the explicit 
knowledge of p(T) for the composite junction. 

The total conductance Gs [normalized to the normal- 
state value Gn — 9192/(91 + 92)] and the Fano fac- 
tor Fs of the S-cavity-N junction are shown in Fig. 2 
as functions of the junction asymmetry and for differ- 
ent couplings between the cavity and leads. For the 
symmetric quantum-point-contact coupling 91/92 = 1, 
the conductance ratio has the minimal value Gs/Gn = 
2(2 - \/2) w 1.17, while the Fano factor is maximal, 
Fs = (v^ + l)/4 PS 0.60. In the highly asymmetric 
limit, Gs/Gn = 2 and Fs = ^gmm/gmax « (Fig. 2, 
solid curves). The vanishing of the shot noise in this 
case is due to the perfect transparency of the dominant 
(the one which is weakly coupled) quantum point con- 
tact. For the case of two tunnel junctions^^''** instead of 
quantum point contacts, the trend is opposite: for the 




0.0 0.2 0.4 0.6 0.8 1.0 



1/ {I + 91/92) 

FIG. 2: Conductance Gs [panel (a)] and the Fano factor Fs 
[panel (b)] of the S-cavity-N junction as a function of the junc- 
tion asymmetry 1/(1 -I- 31/52). Results for three different cou- 
plings of a cavity to the leads are shown for comparison: cou- 
pling by quantum point contacts (solid curves), tunnel junc- 
tions (dashed curves), and coupling by tunnel junction from 
the S side and quantum point contact from the N side (dotted 
curves). Conductance is normalized to the normal-state value 
Gn = 5132/(31 -I- 52), with Qi = {2e^ /h)Ni for quantum point 
contacts and gi — {2e^ /h) J2n J^n^ ^o^' tunnel junctions. Fano 
factors _Fjv for the corresponding normal-state junctions are 
shown in panel (b) for comparison. 

symmetric coupling, Gs/Gn = 1/V2 ~ 0.71 is maximal 
and Fs = 3/4 is minimal, and for the highly asymmetric 
coupling, Gs/Gn = 9rnin/9max « and Fg = 2 (see 
Fig. 2, dashed curves). We point out that these different 
trends can be used to probe the quality of the contacts. 
The dotted curves in Fig. 2 show numerical results for 
the conductance and the Fano factor of a cavity cou- 
pled to a superconductor by a tunnel junction and to a 
normal lead by quantum point contact. We find that the 
transport properties of the system are not affected by the 
type of coupling to the superconductor when the quan- 
tum point contact on the normal side dominates. This 
limit is reached at the conductance ratio 91/92 ^ 5 of 
the tunnel- and quantum-point-contact coupling of the 
cavity to the leads. Therefore, the junction that cor- 
responds to the model of an open chaotic cavity can be 
realized either with two good quality quantum point con- 
tacts from the both sides or it can be an asymmetric 
junction with only the normal-side quantum point con- 
tact of a high quality. The latter is easier to fabricate. 
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and the required asymmetry can be achieved by increas- 
ing the contact area of the cavity to the superconduc- 
tor. Experimentally, the conductance and current noise 
power have been measured recently by Choi et al.,^^ in 
a setup which is very similar to the system we have ana- 
lyzed. As the estimates from Ref. 25 show, it is possible 
to fabricate a high-quality contact between cavity and su- 
perconductor. The measured values of the Fano factors 
Fs = 0.58 ± 0.10 and Fn = 0.25 ± 0.04 across the junc- 
tion in the superconducting- and normal-state regimes, 
respectively, are in agreement with the model of sym- 
metric open chaotic cavity [compare with solid curves in 
Fig. 2 (b)]. However, the measured conductance ratio 
Gs/Gn ~ 0.90 is in discrepancy with the conductance 
ratio Gs/Gn = 1-17 predicted by this simple model. The 
difference may originate from the inelastic quasiparti- 
cle scattering at the disordered superconductor interface, 
nonuniversal correction due to relatively large openings 
of the cavity,^° or dephasing of quasiparticles due to an 
additional lead which is left floating in the experiment. 

At high temperatures ksT ^ \e-V\, current noise 
power is thermally dominated and linear in conductance 
and temperature [see Eq. (24)], as expected from the 
fluctuation-dissipation theorem. Thus, to extract the 
Fano factor from the current noise power measurement, 
it is necessary to be in the low-tcmpcrature, shot noise 
regime fc^T <C \e-V\. Experimentally, this requires high 
bias voltages at which nonlinear I—V characteristics oc- 
cur, especially in a strongly interacting electron systems, 
with a difficulty how to distinguish the shot noise contri- 
bution from the contribution of thermal noise modified 
by nonlinear conductance. 

Finally, we point out the difference between the su- 
perconducting and normal-state asymptotic behavior of 
the third cumulant at high biases. For an asymmetric 
open cavity coupled to normal-metal leads, the slope of 
the voltage-dependent third cumulant is negative, reach- 
ing zero for the symmetric cavity."^* When one lead is 
superconducting, this slope is negative for highly asym- 
metric couplings and positive for symmetric couplings 
[ Fig. 3 (b )], with the crossover at r/ = gmin/9max = 
\/3-|-2\/2- V2 + 2V2 w 0.22. Thus, in the normal-state 
regime the third cumulant changes sign at high enough 
biases, while in the superconducting case this happens 
only if the junction is sufficiently asymmetric. This dif- 
ference originates from the interchannel mixing inside 
the cavity in the presence of the superconducting prox- 
imity effect'^^ and can be attributed to the skewness^^ 
of the Andreev reflection probability distribution func- 
tion pa^Ra) = 2pc{T)dT / cLRa [compare with Eqs. (28)- 
(30)]. For the normal-state symmetric cavity, the trans- 
mission distribution is symmetric - i.e., Pc(^) = Pc{^—T) 
(Fig. 4) - leading to the saturation of the third cumulant 
at high bias. If the jimction is asymmetric, then the 
gap at low transmissions opens at < T < Tq, shifting 
the weight of the distribution towards the open channels, 
Pc{T) < /9c(l - T) for < T < 1/2, and the high-bias 
slope of the third cumulant becomes negative. When 
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FIG. 3: The third cumulant as a function of bias-to- 
temperature ratio [panel (a)], shown for three characteristic 
junction asymmetries: Ci has the maximal positive slope at 

high biases for symmetric coupling rj = Qmin/gmax ~ 1, sat- 
urates at high biases for q ~ 0.22, and has maximal negative 
slope for ?7 « 0.046. The high-bias slope of the third cumulant 
is shown in panel (b) as a function of the junction asymmetry, 
with the normal-state value given for comparison. Effective 
charge is e* = 2e and e* = e for the superconducting and 
normal-state junction, respectively. 



one lead is superconducting, the weight of distribution 
Pa{Ra) for the symmetric cavity is at low values of the 
Andreev reflection probabilities, pa{Ra) > Pa(1 — Ra) 
for < Ra < 1/2 (Fig. 4), leading to the positive high- 
bias slope of the third cumulant. Only for large asymme- 
tries of the junction docs the gap that opens at low Ra 
prevail and the distribution pA shifts towards the open 
Andreev channels and the third cumulant becomes nega- 
tive (Fig. 3). It is interesting to note that the maximally 
negative slopes of the third cumulant in the normal and 
in the Andreev case are approximately equal (if the lat- 
ter are corrected for the effective charge) . From Fig. 4 it 
is seen that indeed the eigenvalue distributions are very 
similar for these values, with the effect of a supercon- 
ductor being the change of a junction asymmetry. We 
believe it would be interesting to confirm these predic- 
tions experimentally. They provide much more detailed 
information on the transmission eigenvalue and Andreev 
reflection eigenvalue distributions, which go beyond the 
information obtained from conductance and noise mea- 
surements. 
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0.0 0.2 0.4 0.6 0.8 1.0 



Ra, T 

FIG. 4: Distribution of the Androov reflection probabilities 
Pa(-Ra) ~ 2pc{T)dT /(IRa (solid curves) and the distribution 
of transmission eigenvalues pc{T) (dash-dotted curves) for an 
open chaotic cavity, shown for different asymmetries of the 
junction r) = gmin/gmax (from top to bottom). 

IV. CONCLUSION 

We have studied the charge transport statistics in co- 
herent two-terminal double junctions within the circuit 
theory of mesoscopic transport. We have shown that the 
system of circuit-theory matrix equations for the Green's 
function of the central cavity can be solved, effectively, as 
a system of scalar equations independently on the type 
of the leads. For an asymmetric cavity coupled to the 
leads by quantum point contacts, the Green's function 
is expressed in a closed analytical form in terms of the 
matrix Green's functions of the leads. The full counting 
statistics and the first three cumulants are obtained for a 
special case of an open cavity between a superconductor 
and a normal metal, at temperatures and bias voltages 
below the superconducting gap. 

The same results can be obtained by applying the cir- 
cuit theory while considering the whole structure as a 
single connector, with the cumulant generating function 
integrated over the distribution of transmission eigenval- 
ues of the composite junction. This approach manifestly 
reveals how the subgap transport in S/N structures is af- 
fected both by the effective charge doubling due to the 
Andreev process and by modification of the transmission 
properties due to electron-hole correlations introduced by 
the superconductor. 

For an open cavity, the Fano factor is enhanced with 
respect to the corresponding normal-state junction, in 
agreement with the recent experimental results by Choi 
et al.,^^ where the high-quality contacts between a cavity 
and superconductor have been made. In comparison to 
the tunnel coupling, the conductance and Fano factor ex- 
hibit opposite trends as a function of the junction asym- 
metry, which can be used to experimentally probe the 
quality of the contacts. The third cumulant is strongly 
affected by the presence of a superconductor. In contrast 
to the normal-state case, in which the third cumulant 



changes the sign at high enough biases, in the case in 
which one lead is superconducting this happens only if 
the junction is sufficiently asymmetric. This difference 
originates from the skewness^^ of the Andreev reflection 
distribution function, which is in favor of closed Andreev 
channels for a moderate asymmetries of the junction. 
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APPENDIX 

Here we show that pi given by Eq. (3) depends only 
on the anticommutator of the Green's functions of the 
leads, pi = pi{{Gi,G2}), under the assumption that it 
is possible to expand Gc in a series over the products 
of G matrices. In this case pi is also a series over G 
matrices, which we denote as pi = pi{G\,G2)- In the 
following we consider only pi, while the reasoning for 
p2 is analogous. In p\{Gi,G2) we separate the term 
ip{GiG2,G2Gi), which contains a sum of products with 
an even number of G matrices, and the term with the 
sum of odd-number products. The latter is of the form 
Gi4>{GiG2, G2G1), where we used G\ = G\ = 1. As a re- 
sult, pi(Gi, 62) = ^(GiG2,G2Gi) + Gi0(GiG2,G2Gi). 
Now we investigate the structure of and (j). First we ex- 
press G1G2 and G2G1 in terms of {Gi, G2J and [Gi, G2], 
and then expand ip = i^({Gi, G2}, [Gi, G2]) in a series 
of [Gi,G2]. Even powers of [Gi,G2] can be expressed 
in terms of {Gi,G2} by using the identity [Gi,G2] = 
{Gi, 62}^ - 4. Thus,^ q> is of the form ^ = d + /3[Gi, G2], 
where a = q;({Gi,G2}) and $ = ${{61,62}) depend 
only on the anticommutator [and hence they commute 
with Gi, G2, and Gc = Gc(Gi,G2)]- The same is true 
for i.e., (f) = d' + /3'[Gi,G2], where a' = Q!'({Gi,G2}) 
and P' = $'{{61,62}). Therefore, 

Pi=a + $[61,62] + a' 61 + /3'Gi[Gi, G2]. (31) 

On the other hand, [Gi,pi] = in accordance with 
Eqs. (3) and (4), and we have 

[Gi,pi] = 2/3Gi[Gi,G2] + 2$' [61, 62] = 0. (32) 

From Eqs. (31) and (32) we find that pi has the following 
structure: pi = a + a' 61. After substituting the expres- 
sion p\6i = a6i +a' back into Eq. (2) from which Gc is 
to be obtained, the term a' vanishes because it commutes 
with Gc- Thus, when solving for Gc we can assume that 
pi = d({Gi,G2}) without loss of generality. A similar 
consideration holds for p2. 
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